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FOREWORD 


The present report is part of broader effort which concerns an. 
investigation of ESSA VII radiation data for use in long-term earth 
energy experiments, performed under Langley Research Center Contract 
No. NAS1-11871 for the National Aeronautics and Space Administration. 
This report Is one of three companion reports which together constitute 
the final report of phases I and II of subject contract. The other 


two reports are entitled; "An Investigation of ESSA VII Radiation 
Data for Use in Long-Term Earth Energy Experiments," published as 
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STEADY -STATE SOLUTION TO THE CONDUCTION PROBLEM 


OF A SPHERICAL BALLOON RADIOMETER 


INTRODUCTION 

House and Sweet (1973) have proposed a satellite system for 
observing the earth's radiation balance employing spherical balloon 
radiometers. The instrumentation consists of three inflatable, 
spherical structures which are exposed to external irradiances of 
the radiation budget. In the steady-state condition of radiative 
equilibrium, the magnitudes of absorbed external irradiances are 
sensed by internal radiometers mounted on the skin of each balloon. 

The temperatures of the radiometers are monitored as a measure of 
the balloons’ internal irradiances (equal to absorbed external 
irradiances) and then telemetered to earth. These temperatures 
will vary depending on the external optical properties of their 
respective balloons, whose properties are selected to effectively 
discriminate the components of the earth's radiation balance. 

In the discussion of measurement principle cited above, two 
assumptions are made concerning conduction of heat: 

• lateral conduction along the balloon skin represents a 
negligible contribution to the overall power density 
balance on an area element of a balloon. 

* the difference between the inside and outside temperatures 
of an area element is negligibly small. 

The effect of lateral conduction becomes small when the balloon 
radius is many times larger than its thickness. In the case of the 

A 

balloons considered above, the radius is about 10 times the thickness. 
Therefore, it is reasonable to neglect lateral heat conduction in the 


first approximation. 
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The assumption that the difference between the inside and 
outside temperatures of an area element is negligibly small implies 
that the skin is quite thin and/or the thermal conductivity is quite 
large. Materials under consideration will lead to balloon thicknesses 
from about 0.5 to 3.0 mils and possibly thicker for those surfaces 
which are painted. Thermal conductivities of materials will vary. 

In the case of aluminum it is large; however, plastics such as mylar 
and teflon, and paints have lower magnitudes of thermal conductivity. 
Small, but significant temperature gradients may be established across 
the skin for the balloon thicknesses and materials currently under 
consideration. Furthermore, the effect of temperature gradients 
through the skin are enhanced by the fourth power of temperature when 
considering radiant power densities according to Stefan-Boltzmann’s 
Law <f T 4 . 

The purpose of this investigation is to assess the conduction 
problem mathematically in order to determine its impact on measurement 
accuracy. It is desirable to develop the problem in terms of the 
magnitude of irradiant sources, balloon thickness and thermal conduc- 
tivity so that the results may be utilized effectively in the 
engineering design of the satellite system under realistic irradiant 
environments . 
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PROBLEM DEVELOIMENT 

Consider the measurement process depicted in figure l. The 

spherical balloon shell absorbs radiant power from irradiant sources 

2 

H over its cross-sectionC), areafT a , assuming the upper limit case 
s 

that all irradiation is parallel. When the balloon is in radiative 

equilibrium with its environment, the absorbed power is balanced by 

2 

an equal loss of power W fi back to space from its surface area 4ir« » 

The magnitude of power loss from an area element depends on the external 
temperature T g according to Stefan-Boltzmann's Law f*T e ^ if the balloon 
is a black body. 

An internal radiometer is located at the inside skin of the 

balloon. This radiometer absorbs radiant power H, from the inside 

- ^ 

of the balloon where T^ . The temperature of this radio- 

meter T r is monitored and eventually telemetered to earth. The point 
in question is whether it is possible to work backwards from an 
observed temperature of the internal radiometer to obtain an accurate 
measurement of irradiant sources? 

Several mathematical expressions hold for the balloon configuration 

under different bounding conduction assumptions. Under conditions of 

infinite conductivity, the balloon will be isothermal. Thus, 

T =T,=T? . Then, 
e i r ’ 


and 


Hs= ^ ^ 

± We- Wi = H-i= 


( 1 ) 


It is evident in equation (1) that H is related to the temperature 

s 

observation T . 

r 

A more general case is to assume that lateral conduction can be 
neglected (i.e., the radius to thickness ratio is very large) and the 
skin thickness s.s quite small. For these conditions temperature gradients 
through the skin are negligible and T^ = T e for an area element. 
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Figure 1. Schematic representation of measurement by a spherical 
balloon radiometer. 
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However, the temperature will vary from one area element to the other. 
Defining the external and Internal average radiant power densities as 

( 2 ) 


and 


V/ e = j|<r7 l+JS 

K-UrV dS 


(3) 


respectively, a similar chain relationship as equation (1) is valid: 

i vE* 131* Hi »«•■&* 

Again the lrradiant sources H are related to the observed tempera- 
ture T r> 

The bounding conduction characteristics assumed in this work are: 

• lateral conduction can be neglected 

• conduction through the skin is a significant part of 
the heat transfer 

When thermal conduction through the skin thickness is slow, T ^ i T g 

and H. ■ W. ^ W by some amount (see Figure 2a). This condition will 
lie 

lead to an error in measuring , utilizing the observation T f . 


Defining an error parameter Q m 2 (W^ - W. ) , the error in the lrradiant 
sources H , defined by £ , is related by 

6 = Z S <6) 


where 


The factor of two in equation (6) is for convenience in later mathe- 
matical development (see equation 9). The goal of the solution to the 
balloon conduction problem is to formulate the error 6 in terms of 
the magnitude of irradiant sources H^, the the thickness of the balloon 
and the thermal conductivity of balloon materials. 
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Figure 2. Schematic representation of (a) balloon conduction problem 
and (b) approximation used in problem solution. Note: the 

irradiance on slab I is 4 H , providing the same average mag- 
nitude of irradiance per unit area as that on hemisphere I. 
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SOLUTION 


The balloon conduction problem can be represented by the diagram 
in figure 2(a). In this case, the sphere is depicted as tvo separate 
hemispheres; one hemisphere (I) is exposed to ir radiant sources H g , and 
the second hemisphere (II) is in the shadow of the first. Since lateral 
conduction is negligible, the transfer of heat from the first to the 
second hemisphere is accomplished solely by radiative transfer within 
the balloon interior. The geometry in figure 2(a) is approximated 
nicely by considering each hemisphere as a semi-infinite slab of similar 
physical properties. This representation is shown in figure 2(b). 

In the slab configuration, irradiation impinges on the outer 
surface of slab I A portion of the radiant power is conducted through 
slab I, radiated to the inner surface of slab II, and, after conduction, 
exits from the outer surface of slab II. The radiometer views both 
slabs equally in figure 2(b) and in a manner similar to viewing both 
hemispheres in figure 2(a). Expressions relating various parameters 
discussed previously are as follows: 

U4 = WL±^ s i.H s 


= Hi* ft* < s > 

S- + u k) (9) 


In figure 2(a), energy (heat) flows per unit area and time from 
the illuminated side, represented by temperature T^, through slabs 
I and II, and is radiated out to space at an effective temperature T^. 
For this energy flow to maintain itself, the relationship between 
temperatures is that T^> T 9 y T^ > T^. 
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A chain relationship depicting this energy flux (power) density is 

t H« ~ = 7 ft- £ )= Hr ^ (10) 

where is the ratio of skin thickness to the thermal conductivity 
of the material in units of m/(watts/(m - K)). The expression 




(Tj - T^) in equation (10), for example, is the power density 


conducted through slab I in units of watts/m . For simplicity the skin 
material is assumed to be homogeneous with an average thermal con- 
ductivity. 

Referring again to figure 2(b), in the condition of radiative 

equilibrium, each of the four surfaces will be at temperatures T^ 

— 4 

and corresponding radiant flux densities ■ <j T j . The thermal 
configuration of the system will depend only upon the magnitude of 
H g , the plate thickness and the thermal conductivity of the materials. 
In this regard, four equations can be written using equation (10) 
to describe the steady-state condition of the system. 

tH s -V, = W s 


(V-TL)=r(kH,-U) 

(T t -z)=y(w,-w t ) 
ui = wt 


(lla) 

(lib) 

(11c) 

(lid) 


Equations lla and lib describe the radiant power balance and the heat 

conduction, respectively, for slab I. Similarly, equations 11c and lid 

describe conditions for slab II. 

The problem at hand is now one of finding the simultaneous 

solution of the system of equations in set (11) with H and as 

s 

adjustable parameters. Since these equations are non-linear in tempera- 
ture, a simple algebraic solution is not possible. One way of approaching 
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T ’ r r 


| r — 1 


this type of problem is to employ a perturbation technique Involving 
a Jacobian operation which leads to an iterative solution to the 
equations. A further simplification results in the realization that 
it is not necessary to find a solution which explicitly determines all 
unknown quantities -- it is necessary only to find the parameter in 
equation (9) . 

The procedure for solution would be as follows: 

1. Choose a ^such that the trivial solution (i.e., y” =0) for Ls 
found. 

2. Perturb the initial choice to the general case (i.e., t> 0). 

3. Using the choice in 1. and the required perturbation as Inputs 

to the Jacobian technique, find a solution for f in terms of 

and H . 
a 

First assume that the material has an infinite conductivity (i.e., 


“ 0). Then the equations in set (11) become 


°=$m= 

V^+Vs 

(12a) 



(12b) 



(12c) 


(12d) 

Under this condition. 



implying that 

%~K 

(13a) 

implying that 

& 

II 

(13b) 

or, Wj . 

Therefore, 



and no error in the measurement results. 
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^ r 


Now perturb the solution above to the realistic case for 
7^ ^0. la doing so equations in set (12) become 

o=^. = fw»- M- + 

o=a*= (V-l)- * (t H s - tf) 
o=3, = (r t -T.)-r (V..-V,) 

0=3* = U 4 - U/ 3 - 


(14a) 

n/.h) 


H4r) 
n 4d) 


The equations in set (14) have been changed by the following amounts 

43 , = a 9* = o= 

A^x = ~n'(z 

&P=-y(W x - W s )=$ s -3.j 

where Aj; , is defined as the required perturbation of the i C *' equation. 
Also for r>o assume that 

0=(Vi <-A ld,)+ ( W*+AK4.)-(U(.+4 t*i)- (ltfj+4 Id*) os) 

Subtracting equation (IS) from equation (9) gives 

S- — [(& W, 4 A U/*) - (A + A W>j] „ 6) 


where the ^W^'s correspond to the changes in the power densities required 
to solve set (14). The ^ W.'s are found using the Jacobian technique. 

The Jacobian operation on set (14) gives, in matrix form, 

J*(W) AW = A<f 

- ty/Wj 


where 


(J<V^ = 
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Incorporating information in equation (10) and set (13), equation (17) 


can be expressed as four simultaneous equations 

+ A*/j = O (l*fl) 

(*+l/+<rT*) &vJ, -(l/t-rT, 3 )Mi=-r(iH s -V/,) ohm 

-r &Wx +(r+ t/+el?)&\ <! s -(i/+r T/)6H^l(Wr tfj) „ *•> 

AK^-AU/j-A^O "««) 


If equations (18a) and (18d) arc added, 

-AW,-4KV=0; Atf,= -AW* 

Therefore, the solution for (^* in equation ( 16 ) reduces to 


<f= (AWi+Att/j) 


(19) 


(20) 


The remaining portion of the problem involves the determination of 

(A W 2 + A W j) from equation set (lb). However, the solution includes variables 

W, and T, in addition to H and Y • It is desirable to 

j j 8 

express the solution only in terras of the adjustable parameters 
and Y • Tl»e following identities follow from equation (10) and 
set (13) which aid in the elimination of and T^ . 

H s = 3 W, ami T l 3 =(H s /3<r) 3/ ‘ 1 ' 

Hs= 6 W/* <•** T* = (Hs/6r) ^ 

The expression for the error parameter Q is 

<f= 0.0303 M?* ( I.Sir 'n/+-hUO.t) (2D 


when a value of £ = 5.6697 x 10 8 watts/nZ-K." 4 is used in the evaluation. 


2 ..A 
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Equation (21) reduces to the compact expression 

<f= 0.0018 * H s 7/lt 


( 22 ) 


3/4 

since the term in the denominator is insignificant for reasonable 

values of H < 2000 watts/m and Y<10 m"-K/watl. From equation (6), 

a 

the error in the observed irradiant sources becomes 


6*0.0036 



(23) 


Equation (23) can also be expressed as a percent error of H . 

• i 8 

Emtotr 100- 036 y /£* (%) (24) 

This concludes the development of a steady-state solution to the 
conduction problem of a spherical balloon radiometer. The impact of 
these observational errors as related to earth radiation balance 


measurements are discussed in the results section below. 


RESULTS 

The solution to the conduction problem in equation (23) indicates, 

as would be expected from the physical characteristics of the problem, 

that the observational error of a balloon radiometer increases as the 

magnitude of irradiant sources increases, as the balloon thickness 

increases and as the thermal conductivity of the skin materials decreases. 

Quantitatively £, increases as the 1.75 power of H and is directly pro- 

s 

portional to the parameter Y . Thus, the observational error increases 

linearly with balloon thickness and inversely with the thermal conductivity 

of the materials. Consequently, skin materials such as alumimm which are 

good thermal conductors have reduced observational errors, whereas, plastics 

which are poor thermal conductors have larger observational errors compared 

to alisninun. Also the solution indicates that thin skinned balloons are a 

desirable engineering design goal. 

Further examination of equation (23) indicates that the errors in 

observations comprise a family of straight, parallel lines on a log-log 

plot of 6 vs. Y* , having slopes of one. Figure 3 presents cuch a plot 

for representative values of H and In order to utilize figure 3, 

s 

all one has to do is compute the ratio of skin thickness to average thermal 

conductivity, and then examine the variation in the observational error for 

different magnitudes of irradiant sources associated with components of the 

earth radiation balance. As an example, assume the parameter ^ - 5 x 10 "* 

in magnitude, the observational error during typical nighttime observations 

2 2 

(H = 240 watts/m ) would be .026 watts/m , and during daytime observations 

s 

2 2 

(H = 1800 watts/m ) is .90 watts/m . 
s 
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Figure 3. Variation in the observational error c with > for different 
magnitudes of irradiant sources H , watt/n.^. 
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In a similar manner, equation (24) can be plotted on a log-log 
graph of 7o error vs. *jT . The graph in figure 4 illustrates such a 
plot for representative values of H . This graph can be utilized in 
the same manner as that in figure 3. The percentage errors associated 
with the examples above for nighttime and daytime observations are ,017a 
and .057>, respectively. 

The balloon radiometers discussed by House and Sweet (1973) consist 
of thin, laminated layers of aluminum and plastic materials which are 
joined by an adhesive. Three balloons are employed in their satellite 
system having different external optical properties. One balloon is 
painted black, the second is white (either a paint or second surface 
mirror) and the third is an aluminum surface. These balloons may be 
denoted by the optical colors as black, white and aluminum. 

Estimates of the observational errors caused by thermal con- 
duction through the balloon skin were made, using the solution developed 
in this work. The thermal properties and construction characteristics 
are needed in order to determine these errors. Table I is a list of 
thermal conductivities for paints and various balloon materials. Table 
2 summarizes the construction details of realistic balloons under con- 
sideration. Weighted average thermal conductivities were computed from 
the coefficients in Table 1, weighted by the thickness of each material. 
Aluminum is neglected in the calculations since the conductivity is large 
compared to other materials and, therefore, contributes little to the 
observational error. Corresponding magnitudes of ”JS* were computed for 
each of the balloon configurations in Table 2. 

It should be. borne in mind when using the graphs in figures 3 and 
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Variation in the percent observational error with y for different 
aagnitudes of irradiant sources H , watt/iL. 
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TABLE 1 


Conductivities of Balloon Materials 


Coefficient of Thermal 
Conductivity k, watt/(m-K) 


0.17 

0.17 

221.0 

0.15 

0.21 

0.20 


0.20 


TABLE 2 


External 
Balloon Color 


Aluminum 


Black 


White 

Painted 


White 

Anodized A1 


White 
2^1 Surface 
Mirror using 
Teflon. 


* 

Aluminum is 
conductivity. 


Thickness of Construction Materials and Conductivity 
Parameters of Typical Balloon Radiometers 


Material 

Thickness 
d ) mil 

Average Conductivity 
k, (watt/ (m-K)) 

Parameter 

"8* , ( (m 2 -K) /watt) 

Aluminum* 

0.5 



Cap tone 

0.2 

0.17 

1,4 x 10 4 

Adhesive 

0.2 



Mylar 

0.5 



Black Paint 

3.0 


-4 

Aluminum* 

0.5 

0.17 

5.5 x 10 

Adhesive 

0.2 



Mylar 

0.5 



White Paint 

5.0 



Aluminum* 

0.5 

0,19 

7.6 x lO- 4 

Adhesive 

0.2 



Mylar 

0.5 



Aluminum* 

2.0 



Adhesive 

0.2 

0.16 

1.1 x 10~ 4 

Mylar 

0.5 


• 

Teflon 

2.0 


_ . _ 4 

Aluminum* 

0.5 

0.20 

3.4 x 10 

Adhesive 

0.2 



Mylar 

0.5 




neglected in the calculations because of its large coefficient of thermal 
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4 that the theory developed herein applies to a perfectly black balloon. 

The balloons listed in Table 2 have varying optical properties both 
externally and internally. The black balloon raoBt closely approaches 
the Idealized situation of the solution to the conduction problem in 
equation (23) . The white balloons will absorb about one fourth as much 
direct and reflected solar radiation as the black balloon. However, 
the thermal emitting properties of the white balloons would be similar to 
that of a black balloon. 

The aluminum balloon is considerably different than either the 
white or black balloons. An aluminum surface absorbs about the same 
portion of solar radiation as the white balloon. But the capacity of an 
aluminum surface to loose heat through thermal emission is about five 
percent that of a black balloon. 

Surfaces exposed to solar radiation frequently are classified by 

their external shortwave absorptivity to longwave emissivity ratio, termed 

the alpha over epsilon ratio. For the surfaces under consideration, the 

black balloon has an optical ratio of about 1.0, the white balloon about 

0.25 and the aluminum balloon about 4.0 in magnitude. As a result, the 

white balloon acquires an equilibrium temperature colder than the black 

balloon, and the aluminum balloon is warmer than the black balloon. 

One approach in the application of the theory to a balloon with 

typical optical properties is to scale the incident irradiance H by the 

s 

corresponding alpha over epsilon ratio. Using the above ratio magnitudes, 

2 

a value of irradiance H equal to 1600 watts /m for the black balloon w^uld 

s 

2 

be increased to 6400 watts/m for the aluminum balloon and reduced to 400 
2 

watts/m for the white balloon. This approach is a first order approximation 


which can be justified for purposes of this preliminary analysis. 


Observational errors for the balloon configurations of figure 2 are 

2 

summarized in figure 3. Night and day values of 240 and 1800 watt/m were 
assumed, respectively, as reasonable magnitudes of irradiance II (see 
black balloon in table 3) . These magnitudes are scaled by the alpha over 
epsilon ratios above and are listed in table 3, for reference. Observa- 
tional errors for the black and aluminum balloons are larger than those 
for the white balloon even though the parameter in table 2 for the 

white balloon is larger than that for the other two balloons. During 

daytime observations, the error magnitude for the white balloon is 0.12 
2 

watt/m , whereas, the errors for the black and aluminum balloons may be 

2 2 

as much as 0.98 watt/m and 2,84 watt/m , respectively. In terms of 

percent errors of H , all values are less than 0.05 percent. 

s 

It should be kept in mind when analyzing the results that the assump- 
tion is made where the incident irradiances all come from the same direction. 
Such an assumption maximizes the observational error. In reality, a 
balloon radiometer in orbit is exposed to irradiances from various di- 
rections and of different magnitudes, depending on the relative positions 
of uhe satellite, earth and sun. When these Irradiances are incident on 
different portions of the balloon skin, the observational error is reduced 
since the surface area of the balloon is at a more uniform temperature. Thus, 
the estimates of observational errors in table 3 are larger than would be 
expected under normal operating conditions in orbit. 

In summary, the results of this analysis indicate that observational 
errors during the day are at least an order of magnitude larger during 
daytime observations than those at night. Observations by all three white 
balloon configurations are acceptable both during daytime and nighttime 




TABLE 3 


Estimated Observational Errors of 
Typical Balloon. Radiometers 


Scaled Irradiances 
H , watt/m^ 


Observational ^ 
Error e, watt/m 


Observational 
Error in (%) 


Night - 240 


Night - 960 


Day - 450 

Night - 60 
Day - 450 

Night - 60 
Day - 450 


0.01 

0.03 

0.001 

• 0.001 

0.02 

0.004 

0.002 

0.003 

0.05 

0.01 
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modes of operation. In the case of the black balloon, observational errors 

2 2 
are less than 1,0 watt/m during the day and probably about 0.5 watt/m for 

realistic orbital conditions. These errors are acceptable for purposes of 

earth energy budget observations. 

2 

Finally, errors for the aluminum balloon are less than 0.1 watt/m 

2 

at night, but may be as large as 2.8 watt/m during daytime operation. 

The latter magnitude may not be acceptable for accurate observations of 
the earth energy budget. However, it should be remembered that the 
solution to the conviction problem developed in this report applies to 
black optical properties, considerably different than an aluminum surface. 
The results here suggest that a more detailed model of the conduction 
problem should be developed to determine a better estimate of the obser- 
vational error for an aluminum balloon radiometer. 
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